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(54) Title: DNA MARKERS FOR MEAT TENDERNESS 

^ (57) Abstract: A method for assessing the tenderness of meat from an animal, comprising the step of testing the animal for the 
presence or absence of a genetic marker selected from the group consisting of: ( 1) an allele of the gene encoding calpastatin (CAST) 

Q associated with peak-force variation or genetic variation located other than in the CAST gene which shows allelic association with 
the CAST allele; and (2) an allele of the gene encoding lysyl oxidase (LOX) associated with intron compression of the semitendinous 
muscle or genetic variation located other than in the LOX gene which shows allelic association with the LOX allele. 
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DNA MARKERS FOR MEAT TENDERNESS 
Technical Field 

The present invention is concerned with genetic 
markers for meat tenderness in animals, and with methods 
5 and oligonucleotide probes for assessing meat tenderness 
in said animals, and a kit for this purpose. The 
invention is useful for the selection of animals which 
show desirable traits in meat tenderness either for 
breeding or to select animals destined to be slaughtered 
10 for food. 

Background Art 

Meat tenderness is an important issue for 
consumers, and one which can influence demand sufficiently 
15 for an especially tender meat to command a premium price 
in the marketplace. The physiological change in muscle 
structure during the postmortem period is complex but 
clearly seems to be at least one factor in meat 
tenderness. The calpain/calpastatin system is an 

2 0 endogenous, calcium-dependent proteinase system, theorised 

to initiate in vivo muscle protein degradation. 
Calpastatin appears to inhibit calpain activity and 
therefore may be assumed to have a role in meat tenderness 
through the regulation of postmortem proteolysis. In 
25 particular, calpain is response for the breakdown of 
myofibril protein, which is closely related to meat 
tenderness . 

It might therefore be suspected that calpastatin 
activity will be related to meat tenderness. Indeed, an 

3 0 increase in postmortem calpastatin activity has been 

correlated to reduced meat tenderness. Nevertheless, 
despite such observations, no clear link between the CAST 
gene, which encodes calpastatin, and meat tenderness has 
been established. 
3 5 For example, Lonergan et al . (1995) undertook a 

restriction fragment length polymorphism analysis at CAST 
and failed to find an association with either calpastatin 
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activity or tenderness in cross bred offspring of sires 
from eight breeds. Chung et al. (1999) measured 
calpastatin activity, Warner-Bratzler Shear Force and 
myofibril fragmentation index in forty-seven purebred 
5 Angus bulls. However, they concluded that "PCR single- 
strand conformation polymorphism analysis of the 
calpastatin gene was not useful for prediction of 
calpastatin activity, myofibril fragmentation index or 
meat tenderness". 

10 It is long known that one of the actions of lysyl 

oxidase (LOX) is to initiate crosslink formation at an 
early stage in collagen f ibrillogenesis (e.g. r Cronlund et 
al., 1985). The action of lysyl oxidase is intensively 
studied with hundreds of publications on a variety of 

15 aspects of its importance in cancer (Giampuzzi et al . , 
2001), the vasculature (Nellaiappan et al.) and other 
tissue and organ systems. 

Variation at the gene itself has not been 
associated with differences in beef tenderness although 

2 0 LOX has always been seen as a strong candidate on 
biochemical grounds for a gene contributing to the 
collagen component of tenderness. Analysis of genetic 
linkage has implicated the genomic region that includes 
LOX in linkage analysis of family variation in adhesion 

2 5 and instron compression of the semitendinosis muscle 

(STADH and STIC; Drinkwater et al., 1999). However, LOX 
itself has not been associated with these measures of 
tenderness through the study of population associations. 

Meat tenderness is a complicated trait because 

3 0 there are many sources of variation that affect postmortem 

meat tenderisation. Some of these are non-genetic effects 
such as the age of the beast, the nature of its feed, 
degree of stress prior to slaughter, carcass chilling, 
postmortem ageing time and cooking and testing methods. 
35 It has been suggested (e.g. Koohmaraie (1994)) that 

approximately 30% of the variation in tenderness in meat 
can be explained by additive gene effects within a single 
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breed , and that approximately 70% of the variation is 
explained by environmental and non-additive gene effects. 
In the Lonergan study the cattle were slaughtered at just 
over 1 year of age (43 0 days), the sample contained only 
5 83 animals of random peak-force values, and the sample 
consisted entirely of crosses between various taurine 
breeds. Likewise, in the Chung study purebred Angus bulls 
only 280 days of age were used. In addition, in neither 
study were the animals selected for extreme peak- force 
10 values, and it therefore seems that environmental and. non- 
fixed genetic effects may have contributed to the failure 
to identify any genetic linkage between the CAST gene and 
meat tenderness. 

15 Summary of the Invention [Revise this once claims are 
settled] 

Through using a protocol designed to reduce the 
influence of fixed genetic and environmental effects, the 
present inventor was unexpectedly able to show allelic 

2 0 association between the CAST and LOX genes and meat 

tenderness. In general terms, therefore, the present 
invention is concerned with genetic markers for meat 
tenderness in animals killed for meat which are 
polymorphisms of the CAST and LOX genes or polymorphisms 
25 which show allelic association therewith. 

Accordingly, in a first aspect of the present 
invention there is provided a method for assessing the 
tenderness of meat from an animal, comprising the step of 
testing the animal for the presence or absence of a 

3 0 genetic marker selected from the group consisting of: 

(1) an allele of the gene encoding calpastatin 
(CAST) associated with peak-force variation or 
genetic variation located other than in the CAST 
gene which shows allelic association with the 

35 CAST allele; and 

(2) an allele of the gene encoding lysyl oxidase 
(LOX) associated with variation in instron 



WO 02/064820 PCT/AU02/001 22 

- 4 - 

compression of the semitendinosis muscle or 
genetic variation located other than in the LOX 
gene which shows allelic association with the 
LOX allele. 

5 Preferably, the allele tested for is located in 

the 3' UTR of CAST, and is typically CAST 3 D/E allele 1, 
having the following partial DNA sequence: 

catttggaaaacgatgcctcacgtgttcttcagtgttctgatttctcat 
gacccctttcctcttGgacttgtgggactgtgtttgatgtttccctgggttgttgttt 
10 ataagtcagtcataaAatactgtgcattgggcacatgtctcctcttgagctgctaatc 
gtaga (SEQ ID NO:l) , 

CAST3 D/E allele 2, having the following partial 
DNA sequence: 

catttggaaaacgatgcctcacgtgttcttcagtgttctgatttctcat 
15 gacccctttcctcttAgacttgtgggactgtgtttgatgtttccctgggttgttgttt 
ataagtcagtcataaAatactgtgcattgggcacatgtctcctcttgagctgctaatc 
gtaga (SEQ ID NO: 2) 

or CAST 3 D/E allele 3, having the following 
partial DNA sequence: 
20 catttggaaaacgatgcctcacgtgttcttcagtgttctgatttctcat 
gacccctttcctcttAgacttgtgggactgtgtttgatgtttccctgggttgttgttt 
ataagtcagtcataaTatactgtgcattgggcacatgtctcctcttgagctgctaatc 
gtaga (SEQ ID NO: 3). 

Reduced meat toughness is selected for by 
25 rejecting animals with the "11" and u 12" genotypes and 

accepting animals with the "22" or "23" genotypes. In the 
sequences given above, the allelic difference is 
highlighted with a capital letter. These three alleles in 
the D/E DNA fragment are due to two SNP (single nucleotide 
30 polymorphisms). The first SNP is at base 2655 of Genbank 
sequence L14450, which is the same as base 2959 of Genbank 
sequence AF159246; it is a G to A change so that allele 1 
has G and alleles 2 and 3 have A. The second SNP is an A 
to T change 58 base pairs 3' to the first SNP. Since only 
35 three alleles have been noted for this region, with 2 
SNPs, it implies that there is complete linkage 
disequilibrium between allele 2 and allele 3, and allele 
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3 may have evolved from allele 2. This is expected since 
they are 58 base pairs apart. For predictive purposes, a 
test of the second SNP which gives a result of allele 3 is 
equivalent to a test of the first SNP giving a result of 
5 allele 2. This is consistent with the peak force values 
of animals that are '23' heterozygotes, all of whom have 
low peak force values. While not wishing to be bound by 
theory, it is believed that these polymorphisms are linked 
to a mutation in or near the calpastatin gene (perhaps in 

10 the promoter or an intron) which results in reduced 
calpastatin expression or activity. 

A further polymorphism has been identified in the 
5' UTR of the CAST gene and other polymorphisms which 
exhibit allelic association with the polymorphism of the 

15 3' UTR, and therefore also act as genetic markers for the 
tenderness characteristics described above, may also be 
present at least within the genomic DNA embracing the 
coding region of the CAST gene and the 5 f UTR and 3' UTR 
regions of that gene. In addition, where there has been a 

20 recent reduction in population size for a species, 

particular haplotypes of individuals will be relatively 
over-represented. If insufficient time has elapsed to 
cause allelic association to decay, there will be linkage 
disequilibrium even for alleles which are far apart. 

25 Livestock species such as cattle have been domesticated 
from a relatively small pool of wild ancestors in recent 
times , and therefore in these species allelic association 
is found between alleles that may be remote physically. 
Thus, it may be expected that regions of genetic variation 

3 0 that are outside the CAST gene will also show allelic 
association with the polymorphisms in the CAST gene 
described above, and therefore will be suitable genetic 
markers for the characteristic of peak-force variation. 
Hence, these polymorphisms may also be used to assess meat 

35 tenderness. 

In particular the CASTS microsatellite 
polymorphism (Nonneman et al, 1999) has been found to be 
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useful as a genetic marker for meat tenderness. Of the 
more common alleles, alleles 7 and 9 have been found to be 
associated with low peak-force and allele 3 to be 
associated with high peak-force. 
5 Therefore, the invention encompasses, in 

preferred embodiments, the further step of testing for the 
presence or absence of one or more additional genetic 
markers such as alleles of the gene encoding calpastatin 
associated with peak-force variation, in particular, with 

10 testing for the presence or absence of CAST 5 allele 7 or 9 
and/or the presence or absence of CAST 5 allele 3. The 
most favorable results when the presence of CAST D/E 
allele 2 has been established is to have CASTS allele 7 or 
allele 9 present also, whereas the benefits of the 

15 presence of CAST D/E allele 2 are to some degree 

counteracted if the animal also possesses CASTS allele 3. 

The IiOX polymorphism has also been shown to be a 
genetic marker for meat tenderness, and allele 1 or allele 
2 may be tested for. Just as for the CAST gene, allelic 

2 0 association may be exhibited to alleles located outside 
the LOX gene. 

According to a second aspect of the present 
invention, there is provided a genetic marker for meat 
tenderness in an animal which is a polymorphic form of the 

2 5 CAST gene, being the CAST3 D/E polymorphism or the LOX 

polymorphism . 

According to a third aspect of the present 
invention there is provided an isolated DNA molecule 
comprising the nucleotide sequence set forth in SEQ ID 
30 NO:l, SEQ ID NO:2 or SEQ ID NO: 3. 

According to a fourth aspect of the present 
invention there is provided an isolated DNA molecule 
consisting of the nucleotide sequence set forth in SEQ ID 
NO:l, SEQ ID NO:2 or SEQ ID NO: 3. 

3 5 According to a fifth aspect of the present invention 

there is provided a method for selecting an animal likely 
to yield meat of improved tenderness, comprising the steps 
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of: 

(1) testing the animal for the presence of an allele 
of the gene encoding calpastatin (CAST) 
associated with low peak-force or genetic 

5 variation located other than in the CAST gene 

which shows allelic association with the CAST 
allele and/or for the presence of an allele of 
the LOX gene associated with the low instron 
compression of the semitendinosis muscle or 
10 genetic variation located other than in the LOX 

gene which shows allelic association with the 
LOX allele; and 

(2) selecting animals which have the CAST and/or LOX 
allele and/or genetic variation in allelic 

15 association therewith. 

Advantageously, in order to assess the tenderness 
of meat from an animal and/or select an animal likely to 
yield meat of improved tenderness testing may comprise the 
steps of: 

20 (1) obtaining a biological sample from the animal; 

(2) extracting DNA from the sample; 

(3) amplifying DNA from the CAST or LOX gene and/or 
from regions of genetic variation which show 
allelic association to polymorphisms of the 

25 relevant one of the CAST or LOX gene; and 

(4) identifying the allele present in the amplified 
DNA. 

Typically the allele identified in step (4) is 
one of CAST3 D/E allele 1, CAST 3 D/E allele 2 and CAST 3 
3 0 D/E allele 3 described above and/or the CASTS alleles 
described above. 

Preferably the biological sample is blood, but 
other biological samples from which DNA can be amplified 
may be used. For example, hair root samples, cheek 
3 5 scrapings, skin samples and the like may be used. 

Typically amplification is performed using the 
polymerase chain reaction (PCR), but other DNA 
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amplification methods such as the ligase chain reaction 
are well known in the art, and may alternatively be used. 

Preferably the alleles are identified by 
polyacrylamide gel electrophoresis techniques such as 
5 SSCP, or by other techniques well known to the person 
skilled in the art such as RFLP analysis. 

In a sixth aspect the invention provides an 
oligonucleotide probe for amplification of a genetic 
marker associated with peak-force variation, said genetic 
10 marker being either an allele of the gene encoding 

calpastatin (CAST) or genetic variation located other than 
in the CAST gene which shows allelic association with said 
allele. 

Typically the probe is selected from the group 
15 consisting of: 

castd 5' cat ttg gaa aac gat gcc tea c 3' 
caste 5' tct acg att age age tea aga gga g 3' 
CAST5U1 5 ' - GT AAAGC C GC AC AAAAC AC AC C C AGG - 3 ' 
CAST5D1 5 ' -GTTTCTGGACCCTCTGGATGAGGAAGCGG- 3 ' . 
20 In view of the designation of the primers as 

CASTD and CASTE, the amplified fragment of the CAST gene 
is referred to from time to time as the CAST D/E fragment 
and the polymorphism as the CAST D/E polymorphism. 

According to a seventh aspect of the present 
25 invention there is provided an oligonucleotide probe for 
amplification of a genetic marker associated with 
variation in instrpn compression of the semitendinosis 
muscle, the genetic marker being either an allele of the 
gene encoding lysyl oxidase (LOX) or genetic variation 
3 0 located other than in the LOX gene which shows allelic 
association with said allele. 

Typically the oligonucleotide probe is an 
oligonucleotide probe selected from the group consisting 
of: 

35 LOX K5: 5' tat cac tga tgt caa acc tg 3' 

LOX K6: 5' act cag gca cca aat age tg 3' 
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According to an eighth aspect of the present 
invention there is provided a kit for use in assessing the 
tenderness of meat from an animal and/or selecting an 
animal likely to yield meat of improved tenderness, 
5 comprising oligonucleotide probes for amplification of at 
least one genetic marker for meat tenderness, said genetic 
marker being either an allele of the gene encoding 
calpastatin (CAST) or genetic variation located other than 
in the CAST gene which shows allelic association with said 

10 allele, or an allele of the LOX gene associated with low 
instron compression of the semitendinosis muscle or 
genetic variation located other than in the LOX gene which 
shows allelic association with the LOX allele, and means 
for amplifying DNA. 

15 The primers used to amplify the DNA are the CASTD 

and CASTE primers and/or the CAST5U1 and CAST5D1 primers 
for amplifying the CAST 5 polymorphism. However, other 
primers able to amplify polymorphisms associated with a 
reduction in toughness in meat are envisaged, whether 

20 these be primers which amplify a polymorphism other than 

the CAST3 D/E polymorphism or CASTS polymorphism, or other 
primers able to amplify the CAST3 D/E fragment of CASTS 
polymorphism. 

The methods of the invention may be used both for 

25 the selection of breeding animals and for the selection of 
unpedigreed animals for entry into feed lots. In the 
latter case, the methods of the invention allow for 
animals with unsuitable pedigrees to be excluded from feed 
lots on the basis that highly tender meat is unlikely to 

3 0 be attained with these animals even after a long feed lot 
holdings. Alternatively, such measurements may allow for 
determination of the optimum time to reach maximum meat 
tenderness. The invention is therefore also concerned 
with animals when selected by the method of the invention, 

3 5 their progeny and the use of both selected animals and 
their progeny for breeding, as well as meat from these 
animals. 



WO 02/064820 



- 10 - 



PCT/AU02/00I22 



The methods of the invention are applicable to 
animals including but not limited to cattle and other 
bovids, including water buffalo and bison, to other 
ungulates, including sheep, goats and deer, and pigs and 
5 chickens. 

Throughout this specification, the words 
"comprise", "comprises" and "comprising" are used in a 
non-exclusive sense, except where, the context requires 
otherwise. 

10 It will be clearly understood that, although a 

number of prior art publications are referred to herein, 
this reference does not constitute an admission that any 
of these documents forms part of the common general 
knowledge in the art, in Australia or in any other 

15 country. 

Brief Description of the Drawings 

Preferred embodiments of the present invention 
will now be described, by way of example only, with 

20 reference to the accompanying drawings, in which: 
Fig 1 is a photograph of a single strand 
conformational polymorphism (SSCP) gel which shows 
genotypes for the CAST 3 D/E polymorphism, from left to 
right, 11, 22, blank, 11, 12, 12, 12, 12, 11, 22. 

25 Fig 2 a & b show the distribution of Warner- 

Bratzler peak- force measurements in the two samples of 169 
and 77 animals respectively. Note that extremes were used 
so there is no middle to the distribution. It does not 
imply that the distribution is bi-modal. Note different 

3 0 scales in the figures . 

Fig 3a & b are a plot of the raw Warner-Bratzler 
peak-force measurements against the CAST genotypes. Note 
the gap in the middle due to the use of extremes of the 
distribution. Note the similarity between the 

35 distributions in the two samples. 

Fig 4 is a boxplot of the residual Warner- 
Bratzler peak-force measurements (XI) for each genotype 
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for the first sample. The median quarter and three- 
quarter percentiles, whiskers and outliers are shown. 

Fig 5 shows the distribution of DNA fragment 
sizes for the CASTS microsatellite. Horizontal axis is 
5 the frequency of each allele and the vertical axis is the 
DNA fragment size. Alleles are labelled in increasing DNA 
fragment size so allele ml in this distribution is less 
than 132 bp, m2 < 136 bp, m3 < 138 bp, m4 < 142 bp, m5 < 
144 bp, m6 < 146 bp, m7 < 148 bp, m8 < 151 bp, m9 < 153 

10 bp, mlO < 155 bp, mil < 157 bp, ml2 < 159 bp, ml3 < 161 
bp. DNA fragments were not found in some of the 2 bp 
bins, e.g., in the less than 134 bp bin, and these are 
either extremely rare or non-existent. 

Fig 6 is a box plot of raw LD peak-force values 

15 along the horizontal axis versus CASTS microsat elite 
allele identity along the vertical axis. The boxes 
contain the median value, represented by the dot, a box 
representing the 25 and 75 percentile and whiskers 
indicating the expected range for the distributions, with 

20 outliers indicated by open circles. Care must be taken in 
interpreting this figure since there are some alleles that 
are rare, such as ml, m2 and ml3 (see Figure 5 for the 
full distribution, so interpretations made on those 
alleles are not particularly informative. Note 

25 particularly that this is a distribution of extreme 

values, so the median value will swing from a low to a 
high value if half the samples are high values. 

Figure 7 is a photograph of a single strand 
conformational polymorphism gel showing the genotypes of 

30 LOX from left to right, 11, 11, 12, 12, 22, 11. 

Figure 8 shows the distribution of instron 
compression measurements for the two samples of 166 and 87 
animals combined. Note that extremes were used so there 
is no middle to the distribution. It does not imply that 

3 5 the distribution is bi -modal. 

Figure 9 shows the distribution of adhesion 
measurements for the two samples of 166 and 87 animals 
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combined. Although the sample was selected for extremes 
of instron compression, this has not been translated into 
a series of extreme adhesion measurements. 

Figure 10 is a plot of STIC versus STADH for the 
5 combined sample. Note the non-uniformity caused by the 
selected STIC values. 

Modes for Performing the Invention 

Example 1 - Identification of CAST3 D/E polymorphism 
10 Cattle were chosen from the DNA Bank of the 

Cattle and Beef Cooperative Research Centre located in 
Brisbane, Australia to have as diverse a genetic and 
phenotypic background as possible. Information stored in 
the CRC Database was used to select animals. Animals of 
15 extremes of peak-force were selected, although animals 

with peak- force measures above 12 were excluded since they 
might have confounded peak- force measurements. In 
essence, the procedure was to select cattle in each 
contemporary group which were of phenotypic extreme 
20 measures, to ensure that no sire was represented by a 
cluster of offspring, that all markets and finishing 
regimes were included in each extreme, so that extremes 
were not biased by being representative of a particular 
market or finishing regime. A total of 169 samples were 

2 5 obtained (Table 1) for the first sample. A second sample 

of 77 animals (Table 6) were analysed as a check to 
determine whether the same allelic association could be 
observed in another sample. 

These DNA samples were genotyped for the CAST 

3 0 (calpastatin) D/E DNA fragment using the primers 

castd 5' cat ttg gaa aac gat gcc tea c 3' 

caste 5 ' tct acg att age age tea aga gga g 3 ' . 

The conditions of the polymerase chain reaction 
(PCR) are an annealing temperature of 60 Celsius, 2.5 mM 
3 5 Magnesium chloride, and reagent mixes obtained from 

Biotech International. The DNA fragments were labelled 
via the incorporation of 32 P dCTP into the fragments during 
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the PCR, and the fragments were visualised by 
autoradiography using X-Ray film exposed overnight at room 
temperature. Alleles were scored in numerical order where 
the fastest migrating allele is number 1. 
5 The genotypes were analysed using generalised 

linear models ( GLM ) following the equation peak-force = 1 
+ genotypes nested within fixed effects + error 
implemented via the S-PLUS software. Fixed effects that 
were considered were breed # finish (Australia, Korea, 

10 Japan), contemporary group (cohort), region (pasture v 

grain, north v south) and the covariate of final weight. 
The genotypes were nested within region and breed since 
pure-bred offspring of taurine sires were not pastured in 
the north. The size of the effects associated with 

15 genotype was estimated by the comparison of variances (eg, 
Anders son-Eklund and Rendel, 1993) . To estimate the size 
of effect associated with genotypic substitution, the same 
model was fitted without the calpastatin genotypes. 
Residuals were extracted and compared to the calpastatin 

20 genotypes. These were analysed using an analysis of 

variance to obtain adjusted means for each genotype Plots 
of raw and residual peak-force values against calpastatin 
genotypes were constructed. 

2 5 Example 2 - Analysis of CAST 3 D/E polymorphism 

There are two common alleles (Figure 1) and at 
least one rare allele for the CAST D/E polymorphism and 
both the common alleles are found in all the breeds, 
although there are clear differences in genotype frequency 

3 0 within the breeds. Zebu breeds have a greater frequency of 

the *11' genotype (Tables 2 and 7) than taurine breeds in 
this sample. 

The raw values (Figures 2a & 2b) were then 
plotted against the CAST genotypes (Figures 3a & 3b) and 
3 5 these associations are sufficiently strong to show visual 
associations between peak-force and genotype. The most 
important genetic effect considered in the literature for 
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CAST, breed or taurine versus zebu, has been carefully 
matched so that there are animals of high and of low peak- 
force from each breed in the sample, and breed is not 
expected to be an explanatory variable here, 
5 The analysis (Table 3) of the CAST genotypes 

shows strong, confirmatory evidence of effects of the CAST 
gene or sequences near the CAST gene on peak- force. The 
analysis shows no effect of breed, but since the sample 
consists of individuals of high and low peak- force for 

10 each breed, this was not unexpected. The size of effect 
associated with this polymorphism is approximately 7.9 
percent of the phenotypic variance estimated as a main 
effect, and the deviance associated with CAST genotype 
nested within breed within region is 121.4 (17 df, P = 

15 0.001894). An un-nested interaction term between breed 
and CAST genotype was calculated for this sample, but is 
was not statistically significant. The GLM of the CAST 
genotypes (Table 4) against the residual peak force 
measurements show a statistically significant level of 

2 0 association similar to that of the CAST genotypes 

considered as a main effect (Table 3) rather than when 
they are nested within region and breed. 

A boxplot of genotypes versus the residual peak 
force measurements (Figure 4) shows clear differences in 

25 distributions and the difference between medians of the 
'11' and the % 22' genotypes are approximately 1.2 kg of 
adjusted peak force. The adjusted means from the analysis 
of variance (Table 5) gives a difference of 1.34 kg of 
peak force between the homozygote genotypes. The overall 

30 standard deviation for the residuals is 1.61. 

The GLM of the confirmatory sample of 77 animals 
showed a statistically significant association between 
CAST genotype and peak force, with the *1' allele 
associated with higher peak forces . When the full model 

35 was calculated, none of the factors were statistically 

significant, possibly as a result of the relatively small 
sample size. Terms in the model were dropped one by one 
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using the reduction in AIC as the criterion. All terms 
except the calpastatin genotypes were dropped (Table 8) in 
this automatic procedure, and these show a deviance of 
17.9 (2 df, P < 0.05) explaining 9.5 percent of the 
5 phenotypic variance. This is similar to the 8.9 percent 
found when the CAST genotypes were compared without other 
factors to peak force in the first sample. 

Discussion 

10 The results presented here indicate that genetic 

variation at the CAST gene is important in explaining 
variation for Warner-Bratzler peak- force measurements 
between individuals irrespective of the breed of origin. 
The sample was chosen to control for the effects of breed 

15 and to spread the sample as widely as possible over 

different sire lines, contemporary groups, feeding and 
finishing regimes; care, was taken to ensure that, as much 
as possible, individuals in either extreme were from each 
breed, contemporary group, feeding and finishing system. 

20 In this way, systematic effects of these factors on peak- 
force were controlled so that the effect of the alleles 
would not be due to inadvertently being carried along by 
other factors affecting peak-force values. Indeed, there 
are statistically significant deviations in peak-force due 

2 5 to allelic substitution at this locus even when there is 

no accounting for the other fixed effects. Inspection of 
the raw data show frequency differences within breeds for 
the different genotypes so that the *1' allele is rarer in 
the extreme with lower peak- force values. 

3 0 This association between the *1 # allele and 

higher peak force measurements is confirmed in a second 
smaller sample of extreme animals. These animals are less 
extreme than those in then first sample, they are the 
left-over extremes, and they clearly show not only that 
3 5 the calpastatin genotypes are important but that in such a 
small sample, other factors known to be important are not 
found to be statistically significant. In a well matched 
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sample such as this it is not of concern, since we attempt 
to remove the effects of the other factors as much as 
possible through the choice of samples to analyse. 

The size of the homozygote substitution is 
5 approximately 1.34 kg of peak force for the LD, equivalent 
to 0.83 of standard deviation. This value is likely to be 
overestimated since the extremes of the distribution were 
used, and a proper estimate will require animals chosen at 
random from the full distribution of peak force. 

10 Nevertheless, this is a useful amount of genetic variance 
associated with a single marker and it is expected that 
this marker would be useful in direct DNA marker tests for 
breeding and feedlot streaming. 

The analysis shows no evidence of a breed by 

15 genotype interaction on peak- force, which means that there 
is no evidence that the allele association is different or 
absent in some breeds. This is interpreted to mean that 
there is no heterogeneity in the breeds for the 
association between calpastatin and peak-force. 

2 0 A positive test for allelic association generally 

means that the causative mutation is close to the DNA 
markers. Associations in other studies have indicated 
that allelic association decays at an extremely rapid rate 
so that DNA markers even relatively close to a 

2 5 quantitative trait locus will find no evidence of 

association (e.g., Coleman et al., 1995; Barendse, 1997). 
This indicates that the causative mutation or mutations 
are extremely close to the CAST D/E polymorphism. 

3 0 Example 3 - Identification of CAST5 microstatellite 

' polymorphism 

To determine whether other polymorphisms in the 
CAST gene are associated with tenderness, both of the 
cattle samples (Tables 1 & 6) were genotyped with the 
3 5 CASTS microsatellite polymorphism (Nonneman et al, 1999) 
which was developed from DNA sequence reported earlier 
(Cong et al . , 1998). 



WO 02/064820 



PCT/AU02/00122 



- 17 - 

The primer sequences to amplify this polymorphism 

are 

CAST5U1 : 5 ' - GT AAAGC C GC AC AAAAC AC AC C C AGG - 3 ' and 
CAST5D1 : 5 ' - GTTTC TGGACC C TC TGGATGAGGAAGC GG - 3 ' 
5 with amplification fragments in the range 13 0 - 159 bp, 

sizes determined on an ABI 373 DNA sequencer. Alleles and 
genotypes were assigned based on these size fragments 
leading to 13 alleles and the distribution of allele sizes 
is shown in Figure 5. 
10 Two different sets of analyses were performed. 

In the first, the genotypes at the CAST3 D/E polymorphism 
were compared to the CASTS microsatellite to determine 
whether there were significant associations between the 
genotypes, as a consequence of haplotypes existing along 
15 the DNA sequence. If CASTS ^nd CAST 3 show significant 
haplotypes, since they are on either side of the CAST 
coding sequence, then all polymorphisms for the CAST 
coding sequence are expected to be in linkage 
disequilibrium with either or both of these DNA markers. 
20 In the second, the CASTS microsatellite alleles were 

compared to the LD peak- force measurements to determine 
whether there was any association between CASTS and 
tenderness. 

Haplotypes between CASTS and CAST 3 

Since genotypes of parents of these animals were 
not available haplotypes where determined by analysing 
animals in which one or both of CASTS and CAST 3 had 
homozygous genotypes. The frequency of these haplotypes 
were tabulated (Table 9) . These frequencies were tested 
for heterogenity using a generalised linear model and 
found to be highly heterogenous (Table 10) . This means 
that each allele at CAST 3 D/E is preferentially associated 
with specific alleles at the CASTS microsatellite. 

Association between tenderness and CASTS 

Since CASTS has 13 alleles and hence there are 91 
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possible genotypes, not all of these genotypes will be 
seen in a sample of 240 samples, as in this study, so the 
association was estimated on the alleles. As for the 
CAST 3 D/E DNA marker, the polymorphism was compared to 
5 the raw LD peak-force values (Table 11 a), was examined 
for differences in interactions between breeds (Table 11 
b), and was compared to the LD peak-force values after 
market, cohort, breed and finish effects were accounted 
for (Table 11 c) . In the last of these analyses, CAST5 

10 alleles are nested within finish and breed, as in the 
analysis of CAST3 . 

These analyses show that there is no interaction 
between CAST5 allele frequency and breed on LD peak- force, 
that the association between CAST5 and the raw LD peak- 

15 force values is statistically significant at the threshold 
P < 0.01, but when the CASTS alleles are nested within 
breed and finish, the association has a deviation which is 
0.1 > P > 0.05. The lack of interaction between CASTS and 
breed in explaining LD peak- force means that any 

2 0 differences in gene frequencies between breeds are not 

responsible for the association between CASTS and LD peak- 
force* The association between CASTS and LD peak- force in 
sections a and b of Table 11 indicate that there is some 
evidence for CASTS associated with LD peak-force, but a 

2 5 bias might still exist, which is why the factors such as 

market, cohort, breed and finish are also corrected for. 
Once those factors are corrected, there is a lack of 
strength in the association. In the CAST3 D/E analysis, 
correcting the additional factors improved the evidence 

3 0 for the association, and since the same samples are used, 

we know in which direction the deviations should go. Thus 
the lack of strength probably means that the large number 
of alleles nested within breed and finish, has failed to 
find an association due to the creation of a large number 
3 5 of categories. Larger numbers of alleles are expected to 
reduce the strength of associations purely due to the 
number of categories (cf Terwilliger, 1995). 
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The CASTS polymorphism can be used in conjunction 
with the CAST 3 D/E polymorphism to predict LD peak-force. 
For CAST 3 D/E the ell genotype is associated with higher 
peak force values, the cl2 genotype is intermediate and 
5 the c22 genotype has the lower peak force values. 

Secondly, there is linkage disecjuillibrium between CAST 3 
D/E and CASTS. By examining the table of haplotypes, 
looking at the common microsatellite alleles, CAST3 D/E al 
(allele 1) is most often associated with CASTS m3 (allele 

10 3) with low abundances for m7 and m9 . On the contrary, 
CAST 3 D/E a2 (allele 2) is most often associated with 
CASTS m9, with a similar large association to m3 and a 
lesser but 1 still significant association with m7 . 
Inspection of Figure 6, a plot of raw LD peak-force values 

15 for each CASTS microsatellite allele, indicates that CASTS 
m7 and m9 have lower peak force values while CASTS m3 has 
higher peak force values. Since most of the m3 alleles 
are actually associated with CAST 3 D/E a2 and not al (108 
versus 14), this higher value is not likely to be the 

20 effect of CAST D/E al. Rather it provides a tool to 

refine the assignment of animals to groups, since animals 
selected for having CAST3 D/E a2, so that they might have 
lower peak force values, might still have higher peak 
force values if they possessed CASTS m3 . They are expected 

25 to have a greater likelihood of having lower peak force 
values if they possessed both CAST 3 D/E a2 as well as 
CASTS m7 or m9 . 

Example 4 

30 This example shows the testing of a DNA marker in the 

LOX gene for population associations to STIC and STADH. 
Repeated statistically significant positive associations 
were found between genotypes and both STIC and STADH. 
These indicate that, unusually, the heterozygote may be 

35 one of the extreme genotypes suggesting some 

overdominance. These associations are found in a study of 
6 breeds of cattle with a structure to detect linkage 
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disequilibrium and would indicate that the gene LOX either 
contained or was located near to the genetic factor 
associated with connective tissue strength. 

Materials and Methods 

5 Cattle were chosen from the CRC DNA Bank to have 

as diverse a genetic and phenotypic background as 
possible. Two groups of animals were chosen, the first 
and larger set to test for associations and the second 
smaller set to confirm the polarity of the associations 

10 (cf. Barendse 1997; Barendse et al., 2000). Information 
stored in the CRC Database was used to select animals. 
Animals of extremes of instron compression in the 
semitendinosis muscle were selected. Adhesion measures 
for these animals were also extracted from the database. 

15 In essence, the procedure was to select cattle in each 
cohort which were of phenotypic extreme measures, to 
ensure that no sire was represented by a cluster of 
offspring, that all markets and finishing regimes were 
included in each extreme, so that extremes were not biased 

2 0 by being representative of a particular market or 

finishing regime. A total of 253 individuals were 
selected comprising a first sample of 166 animals and a 
second sample of 87 animals (Table 11) . 

The DNA was genotyped for the LOX (Lysyl Oxidase) 
25 DNA fragment using the primers LOX K5 : 5' tat cac tga tgt 
caa acc tg 3' and LOX K6: 5' act cag gca cca aat age tg 
3'. The conditions of the polymerase chain reaction (PCR) 
are an annealing temperature of 60 Celsius, 2.5 mM 
Magnesium chloride, and reagent mixes obtained from 

3 0 Biotech International. The DNA fragments were labelled 

via the incorporation of 32 P dCTP into the fragments during 
the PCR. The fragments were digested with Hinfl overnight 
at 37 Celsius before separation on gels. The fragments 
were visualised via autoradiography to X-Ray film 
3 5 overnight at room temperature. Alleles were scored in 



WO 02/064820 



PCT/AU02/00122 



- 21 - 

numerical order where the fastest migrating allele is 
number 1. 

The genotypes were analysed using generalised 
linear models (GLM) following the equation STIC = 1 + 
5 genotypes nested within fixed effects + error implemented 
via the S-PLUS software. The same model is used for 
STADH. Fixed effects that were considered were breed, 
finish (Domestic, Korea, Japan), cohort, region (pasture v 
grain, north v south) and the covariate of age. Age was 
10 included since LOX affects cross-linking of collagen and 
cross-linking is an age related process, with cross- 
linking increasing over time* The genotypes were nested 
within region and breed since pure-bred offspring of 
taurine sires were not pastured in the north. The size of 
15 the effects associated with genotype was estimated by the 
comparison of variances (eg, Andersson-Eklund and Rendel, 
1993). To estimate the size of effect associated with 
genotypic substitution, the same model was fitted without 
the LOX genotypes. Residuals were extracted and compared 
2 0 to the LOX genotypes. These were analysed using an 

analysis of variance to obtain adjusted means for each 
genotype . 

Results 

There are two alleles (Figure 7) for the LOX 
polymorphism and both these alleles are found in all the 
breeds, although there are clear differences in genotype 
frequency within the breeds. There is no consistent 
difference between zebu and taurine breeds in frequency of 
the genotypes (Table 12). The Hereford breed differs 
radically in genotype frequencies to all the other breeds 
in the sample. It has high frequencies of genotype y 22' 
while all other breeds have high frequencies of genotype 

The STIC and STADH values are correlated with 
R=0.52 (Figures 8 - 10). The plots indicate that while 
the STIC values show two clear extremes the STADH values 
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have only a long tail and do not show two discrete 
distributions. This reflects that the sample was selected 
only on STIC. 

The analyses (Tables 13 and 14) of LOX against 
5 STIC and STADH show consistent statistically significant 
associations. The first and the second samples as well as 
the combined samples of both STIC and STADH show 
associations to LOX genotypes at P<0.05. The association 
of LOX appears stronger to STADH than to STIC. The 

10 association in the second sample of STADH phenotypes has 
extremely high statistical significance but this may be 
due to sampling in small populations and the cong}Mence of 
extreme phenotypes with particular genotypes. The 
combined analysis of STADH is less extreme than the second 

15 sample but shows confirmatory linkage to LOX (P<0.01). 

Nevertheless, it is clear that these are not 
large genotypic substitution effects and some analyses do 
not show statistical significance. When the LOX genotypes 
were compared to residual STIC and STADH, none of these 

20 associations was statistically significant, whether by 

sample or data combined (Table 15) . Interestingly, some 
of the comparisons show that the heterozygotes are of 
extreme phenotype, opening the possibility of 
overdominance at this locus. This will need to be 

25 confirmed using other polymorphisms at the LOX gene that 
show larger genotypic substitution effects. 

Discussion 

Consistent with those earlier analyses, the STADH 
30 values show greater association to the DNA marker than the 
STIC values, even though the samples are extreme for STIC, 
with STADH values only more dispersed than normal due to 
the correlation between traits (Figure 9). STIC was 
chosen upon which to select extremes rather than STADH. 
35 However, both of these measurements evaluate aspects of 
connective tissue strength, the adhesion measures the 
force required, in crude terms, to pull a muscle apart, 
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the force applied perpendicular to the fibre bundles, 
while the instron compression measures how much the muscle 
can be flattened without being torn or cut. These are not 
perfectly correlated as can be seen by inspection of the 
5 distribution of STIC and STADH values (Figure 10). 

Example 5 

The association between the marker and STIC was 
examined in Example 4 using two batches of extreme 

10 animals. The results show that there are significant 

associations between the genotypes of the marker and STIC 
(instrom compression, P < 0.05)/ and STADH (adhesion, P 
<0.01). The results suggest that the gene LOX either 
contains or is located near the genetic factor associated 

15 with connective tissue strength. 

Because this study was carried out on a relative 
small population (253) with extreme animals only, the same 
marker was tested on different populations to see if the 
association is still valid. 

20 

Materials and Methods 

In addition to the population, there are two 
other groups containing animals chosen from the two tails 
of instron compression (LDIC, 136) and peak force (LDPF, 
25 131) for the LOX gene study. These three extreme groups 
together with 559 non-extreme individuals form the base 
for these analyses on the LOX marker.. A total of 917 
individuals were used for the study (Table 16) . 

Due to the nature of the populations, the 
30 analyses were carried out to the three datasets. 

Extreme animals only (389) . The extreme animals 

from LDIC, LDPF and STPF were pooled together. 

Non- extreme animals (559) . 

Combined data (917). The combination of 1 and 2. 
35 Beside the traits STIC, LDPF and LDIC, a range of 

other traits was also evaluated to see if there is any 
effect of LOX gene on other meat quality traits (Table 
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17). The intramuscular fat measurements from LD__FAT% and 
NIR_FAT% were combined to make a single trait . 

The mixed model procedure (MLX) in SAS (version 
8.0) was used to run the statistical analyses. The fixed 
5 . effects in the model include finish group and LOX marker. 
Sire and contemporary groups are treated as random 
effects. All these effects were nested within individual 
breeds. The statistical model used is as follow: 

Trait = mean + sire within breed + contemporary 
10 group within breed + finish within breed +LOX within breed 
+Carcass weight. 

Contemporary group was defined as the 
combination of herd of origin, cohort and kill code. The 
individuals without electrical stimulation were removed 
15 from the analysis data. Carcass weight is being used as a 
covariate to adjust for the age difference. 

A full contrast model would be performed if a 
significant marker-trait association was identified from a 
mixed model (or GLM) analysis. The purpose of conducting 

2 0 such the test is to further examine the possibility of 

additive or dominance or overdominance effect among the 
genotypes of the LOX marker. The full contrast of 3 
genotypes (11, 12 and 22) is set up in SAS as follow: 

contrast 'Additive Test' lox(bcode) 1 0 -1; 
25 Contrast 'Homozygotel vs Heterozygote' 1-10 

Contrast 'Heterozygote vs Homozygote2' 0 1-1 

contrast 'Dominance Test' lox(bcode) -1 2 -1; 

contrast 'Recessive Test' lox(bcode) -1 -1 2; 

contrast ' Over Dominance Test' lox(bcode) 2 -1 -1; 

3 0 For the extreme population in which the animals 

with extreme phenotypes were genotyped, multi-trait 
logistic regression method (Henshall and Goddard, 1999) 
was applied to take the potential correlation of traits 
into account. The method is regression based,, but instead 
3 5 of regressing phenotype on genotype, the regression is 
genotype on phenotype. This replaces the assumption that 
phenotypes are unselected with the assumption that there 
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was no selection based on genotypes (Henshall and Goddard, 
1998). Prior to using logistic regression method, MLX 
model was used to all data (917 animals) to derive 
predicted values of individual animals. The predicted 
5 phenotype values for the extreme animals after adjusting 
for significant fixed effects were then used for Logistic 
regression analyses. The analyses started with single 
trait logistic regression method and then proceeded to 
multi-trait logistic regression method. 

10 The genotype frequency distribution of the 

marker in different populations is shown in Table 18. From 
the table, it can be seen that the Hereford breed differs 
remarkably in genotype frequencies to all the other breeds 
in the populations. In order to investigate the potential 

15 effect of skewed genotypes of Hereford breed on the 
overall results, a set of additional analyses were also 
pursued to the populations by removing the Hereford 
individuals from the data sets. 

Results and Discussions 
20 Part I. Extreme Animals (Table 19) 

Extreme animals for LDIC 

The first test was conducted to the sample 
containing the selected animals for LDIC (13 6) . The 
25 results from the analysis of variance reveal that there 
was no close association between any genotypes of LOX 
marker and LDIC. The same conclusion was held to other 
meat quality traits. 

3 0 Extreme animals for LDPF 

Like LDIC sample, there was no significant 
variation detected between the LOX marker and any meat 
quality trait in the batch animals selected for LDPF. The 
results are not surprising as the initial QTL for 

3 5 tenderness in CBX experiment was identified in instron 
compression measurement of Semitendinosus muscles. 
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Combined extreme animal data 

Analysis of Variance. As sire effect was 
confounded with other effects , it had to be removed from 
the model and GLM (generalised linear model) was 
5 performed. In this case, contemporary group was treated as 
a fixed effect rather than a random effect. The analysis 
of variance has shown that out of 21 meat quality traits 
tested, STIC and LDL had significant results (P < 0.05). 

Full Contrast Model. The results from full 

10 contrast model are given below. For STIC, it can be seen 
that there was no additive effect between the two 
homozygous genotypes (11 and 22). However, the highly 
significant difference between the phenotypes of 11 and 12 
obviously contributed to the detection of dominant and 

15 overdominant effects. Nothing was remarkable for LDL. 

Logistic Regression. After adjusting for the 
significant fixed effects on all data, logistic regression 
was applied to STIC. Multi-trait logistic regression model 
was also fitted to take the potential correlations between 

2 0 ST measurements into account (STIC, STPF and STADH) . The 
results confirm the findings from the other methods. That 
is, LOX genotypes did have a correlation with STIC. The 
regression co-efficiency between lox marker and STIC is 
shown in the output of logistic procedure (below) . The 

2 5 allele substitution effect of the lox marker could be 

derived from the formulae suggested by Henshall and 
Goddard (1999) based on the total variance of whole data. 
The multi-trait logistic regression test on STIC, STPF and 
STADH has shown that both STIC and STPF had significant 

3 0 effects on LOX gene marker. STPF was marginally non- 

significant in GLM analysis. (Table 20) 

Part II. CRC Non-Extreme Animals 

The non-extreme animals (559) were genotyped 
35 against LOX marker in CRC I and but were not part of the 
animals involved in marker evaluation Phases III. The 
mixed model analyses of variance show that beside STIC, 
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the significant marker-trait association was also detected 
to the intramuscular fat (FAT) and LDPH. However, full 
contrast test to STIC and FAT has failed to pinpoint the 
genotype causing the significant results. In the case of 
5 LDPH, it seems that 22 genotype had an important role in 
determining the outcomes. (Table 21) 

Part III. Combined Data 

When extreme and non-extreme animals were pooled 

10 together, the results from mixed model analysis of 

variance show that again the lox marker was associated 
with STIC (P < 0.05). The significant results were also 
found in STL, which is the measurement of darkness of 
carcass meat colour. However in both cases, full contrast 

15 model had not be able to identify the significant genotype 
sources. (Table 22) 

Part IV. Removing Hereford individuals from the combined 
population 

2 0 In order to test the possible effect of skewed 

distribution of lox genotypes of Hereford breed, the 
additional analyses were also performed to the combined 
data with the removal of Hereford breed. The results 
indicate that the removal of Hereford animals has changed 

2 5 little to the overall significant results of STIC in the 

combined population. From the genotype frequency 
distribution table, it can be seen that the majority of 
Hereford individuals were from the three extreme 
populations except one animal from non-extreme CRC 
30 population. (Table 23) 

The overall results from the investigation of 
LOX gene effect on meat quality traits have been 
consistent across three populations (extreme, non-extreme 
and combined) . That is, there is a strong association of 

3 5 LOX gene marker with the instron compression measurement 

of Semitendinosus muscles (P < 0.05). The significant 
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results from other meat quality traits vary from one 
population to another. 

Industrial Applicability 
5 The invention is useful in allowing selection 

and breeding of animals which yield more tender meat. 
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Table 1 

Characteristics of the first Cattle Sample 

5 Total: 169 





83 


high peak force 




86 


low peak force 


Breeds : 


29 


Santa Gertrudis 




25 


Hereford 




26 


Angus 




27 


Belmont Red 




31 


Brahman 




31 


Shorthorn 


Regions : 


38 


Pasture South 




28 


Pasture North 




57 


Grain South 




41 


Grain North 


Markets : 


72 


Korean 




67 


Domestic 




25 


Japanese 


Cohorts : 


27 


Cohorts 



Median: 5 steers per cohort 
bottom quartile: 2 steers per cohort 
top quartile: 9 steers per cohort 
25 Sires: 112 sires 

Median: 1 steer per sire 

bottom quartile: 1 steer per sire 

top quartile: 2 steers per sire 
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Table 2 

Distribution of CAST genotypes in the breeds in the first 
sample . 

5 



11 12 22 23 

10 Angus 0 7 19 0 

Belmont Red 0 8 19 0 

Brahman 6 13 10 2 

Hereford 0 5 17 0 

Santa Gertrudis 3 5 19 0 . 

-15 Shorthorn 0 4 23 0 
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Table 3 

Associations between calpastatin genotypes (castg) and 
tenderness . 

5 

A. Calpastatin by itself 
Analysis of Deviance Table 

10 

Gaussian model 

Response: peakforce 

15 Terms added sequentially (first to last) 

Df Deviance Resid. Df Resid. Dev F Value Pr(F) 

NULL 155 864.6307 

20 castg 3 70.89899 152 793.7317-4.52573 0.004536025' 



25 

B. Breed x Calpastatin Interactions 



Analysis of Deviance Table 

30 

Gaussian model 
Response: peakforce 
35 Terms added sequentially (first to last) 



40 





Df 


Dev 


iance 


Resid. Df 


Resid. 


, Dev 


F 


Value 




Pr (F) 


NULL 








155 


864. 


, 6307 










finish 


2 


101 


.8455 


153 


762. 


7852 


15 


.97637 


0. 


. 0000008 


cohort 


25 


273 


. 5882 


128 


489. 


.1971 


3 


.43339 


0. 


.0000041 


region 


3 


59 


.6620 


125 


429. 


.5350 


6 


.23940 


0. 


.0005889 


breed- 


4 


10 


.8711 


121 


418. 


. 6639 


0 


.85267 


0, 


.4948672 


castg 


3 


28 


.2709 


118 


390. 


,3930 


2 


.95654 


0. 


.0355308 


breed:castg 


6 


33 


. 4066 


112 


356. 


.9864 


1 


.74682 


0, 


.1166518 



45 
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Table 3 (continued) 

C. Calpastatin genotypes nested with breed and region • 

Analysis of Deviance Table 
Gaussian model 
Response: peakforce 

Terms added sequentially (first to last) 

15 Df Deviance Resid. Df Resid. Dev F Value 

Pr{F) 

NULL 155 864.6307 

finish 2 101.8455 153 762.7852 18.28842 

0.0000002 

20 cohort 25 273.5882 128 489.1971 3.93027 

0.0000005 

region 3 59.6620 125 429.5350 7.14235 

0.0002128 

breed in region 7 26.8943 118 402.6408 1.37983 

25 0.2219802 

castg in (region/breed) 17 121.4139 101 281.2269 2.56498 

0.0018938 



30 
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Table 4 

Analysis of CAST against residual peakforce measurements (XI) 

Analysis of Deviance Table 

Gaussian model 

Response: XI 

Terms added sequentially (first to last) 

Df Deviance Resid. Df Resid. Dev F Value Pr(F) 

15 NULL 154 402.3160 

castg 3 21.97947 151 380.3365 2.90874 0.03654659 



10 



20 



Call: 

glm(formula = XI castg, data = calppftest, na. action = na. omit) 
Coefficients : 
25 (Intercept) castgl castg2 castg3 

0.07693593 -0.4095948 -0.3102917 -0.3504961 

30 Degrees of Freedom: 155 Total;. 151 Residual 



35 



55 



Residual Deviance: 380.3365 



Model from which residuals were calculated 

Analysis of Deviance Table 
40 Gaussian "model 

Response: peakforce 

Terms added sequentially (first to last) 

45 

Df Deviance Resid. Df Resid. Dev F Value Pr(F) 

NULL 160 906.5175 

finish 2 106.0217 158 800.4958 15.52293 0.0000010 

50 cohort 26 290.8387 132 509.6571 3.27558 0.0000057 

finlwt 1 8.0662 131 501.5908 2.36200 0.1269337 

region 3 63.9431 128 437.6478 6.24139 0.0005629 

breed in region 7 24.4323 ' 121 413.2154 1.02206 0.4192678 



glm( formula: peakforce = finish + cohort + finlwt + region/breed, 
data = calppf, na. action = na.omit) 
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Table 5 

Analysis of Variance tables between CAST genotypes and residual 
peak force measures (XI) along with the table of adjusted means 
associated with each genotype. 

Analysis of Variance Table 

Response: XI 

Terms added sequentially (first to last) 

Df Sum of Sq Mean Sq F Value Pr(F) 

15 castg 3 21.9795 7.326490 2.90874 0.03654659 

Residuals 151 380.3365 2.518785 



10 



20 



30 



Tables of adjusted means 



25 Grand mean 

0.076936 
se 0.318703 



castg 



ell cl2 c22 c23 

. 1.1473 0.3281 -0.1932 -0.9746 
35 se 0.5290 0.2479 0.1564 1.1222 
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Table 6 

Characteristics of the second sample of 77 animals. 



10 



15 



20 



25 



Total: 77 

39 high peak force 

38 low peak force 
Breeds: 11 Belmont Red 

11 Hereford 
13 Brahman 

13 Shorthorn 

14 Santa Gertrudis 

15 Angus 
Regions: 24 Pasture South 

12 Pasture North 
21 Grain South 
20 Grain North 

Markets: 35 Korean 

2 5 Domestic 
17 Japanese 

Cohorts: 22 Cohorts 

Median: 3 steers per cohort 

bottom quartile: 2 steers per cohort 

top quartile: 5 steers per cohort 

Sires: 64 sires 

Median: 1 animal per sire 

bottom quartile: 1 animal per sire 

top quartile: .1 animal per sire 
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Table 7 

Distribution of CAST genotypes in the second sample 



5 



Breed Genotype 



11 12 22 

Angus 0 3 12 

Belmont Red 0 3 8 

10 Brahman 3 7 3 

Hereford 0 3 8 

Santa Gertrudis 14 9 

Shorthorn 0 0 13 
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10 



Table 8 

Associations between calpastatin genotypes and LD peak force in 
the second sample. 

Analysis of Deviance Table 

Gaussian model 

Response : ldpeakf orce 

Terms added sequentially (first to last) 

Df Deviance Resid. Df Resid. Dev F Value Pr(F) 

15 . NULL 76 205.9332 

castg 2 17.90313 74 188.0300 3.522925 0.03455689 

Coefficients: 
(Intercept) castgl castg2 
20 5.227591 -0.1205 -0.3719088 

Degrees of Freedom: 77 Total; 74 Residual 
2 5 Residual Deviance: 188.03 



3 0 Single term deletions 
Model : 

ldpeakforce = lslortwait + buttemp + finish + cohort + region + 



3 5 breed + castg 



Final Call: 



40 glm( formula .= ldpeakforce castg, data = calppfr, na. action 
na.omit) 



Coefficients : 
45 (Intercept) castgl castg2 

5.195064 -0.07055556 -0.37438 

Degrees of Freedom: 67 Total; 64 Residual 

50 Residual Deviance: 161.5878 
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10 



15 



20 



25 



30 



35 
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Table 9 

The amount of each haplotype found between the alleles of the 
casts microsatellite and the cast3 D/E SNP on both cattle samples. 

Twenty- six haplotypes were found in animals that are homozygous 
for one or the other locus- 

allele 





cast 3 


cast 5 


amount 




ax 


mx 




2 


ax 




u 




ax 


mj 


14 


A 
*& 


ax 


— >* 
m4 


0 


c 


ax 


m5 


0 


0 


al 


m6 


3 


/ 


ax 


m7 


6 


o 
o 


al 


m8 


2 


9 


al 


m9 


6 


10 


al 


mlO 


5 


11 


al 


mil 


0 


12 


al 


ml 2 


1 


13 


al 


ml 3 


0 


14 


a2 


ml 


0 


15 


a2 


m2 


1 


16 


a2 


m3 


108 


17 


a2 


m4 


1 


18 


a2 


m5 


4 


19 


a2 


m6 


3 


20 


a2 


m7 


42 


21 


a2 


m8 


2 


22 


a2 


m9 


110 


23 


a2 


mlO 


17 


24 


a2 


mil 


6 


25 


a2 


ml 2 


1 


26 


a2 


ml3 


1 
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Table 10 

A heterogeneity test for associations between alleles at CAST5 
with alleles at CAST 3 D/E. 

5 Analysis of Deviance Table 

Poisson model 



Response: score 

10 

Terms added sequentially (first to last) 

Df Deviance Resid. Df Resid. Dev Pr(Chi) 

NULL 25 926.9836 

cast3 1 220.4994 24 706.4842 0.0000000000 

15 casts 12 671.7497 12 34.7345 0.0000000000 

cast3:cast5 12 34.7344 0 0.0001 0.0005161421 



Table 11 

20 Tests for association between CAST 5 microsatellite and LD peak 

force measurements in both cattle samples. 

Part A. Calpastatin by itself 

2 5 Analysis of Deviance Table ' 

Gaussian model 



Response: ldpf 

Terms added sequentially (first to last) 

Df Deviance Resid. Df Resid. Dev F Value Pr(F) 
NULL 491 2266.640 

malll 12 136.0104 479 2130.629 2.548113 0.002843108 



Part B. Breed by calpastatin interactions 



Analysis of Deviance Table 
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10 



15 



20 



Table 10 (Continued) 



Gaussian model 



5 Response: Id 

Terms added sequentially (first to last) 





Df 


Deviance 


Resid. 


Df Resid. 


Dev 


F 


Value 




Pr(F) 


NULL 






491 


226 6 


.640 










market 


2 


218.4756 


489 


2048 


.164 


42. 


35776 


0 


.0000000 


cohort 


26 


706.3980 


463 


1341 


.766 


10. 


53504 


0 


.0000000 


finish 


3 


97.7159 


460 


1244 


.050 


12. 


63003 


0 


.0000001 


breed 


4 


23.5779 


456 


1220 


.472 


2. 


28562 


0 


.0594827 


casts 


12 


65.7617 


444 


1154 


.711 


2. 


12497 


0 


.0146070 


I: cast5 


26 


76.7170 


418 

J 


1077 


.994 


1. 


14414 


0 


.2865614 



Part C. Calpastatin genotypes nested with breed and finish (region) 
Analysis of Deviance Table 



Gaussian model 



Response: Id 



25 Terms added sequentially (first to last) 

Df Deviance Resid. Df Resid. Dev F Value Pr(F) 
NULL 491 2266.640 

market 2 218.4756 489 2048.164 42.90121 0.00000000 

cohort 26 706.3980 463 1341.766 10.67020 0.00000000 

30 breed 4 30.7280 459 1311.038 3.01697 0.01799825 

finish %in% reed 6 111.2067 453 1199.832 7.27908 0.00000022 

casts %xtl% (breed/finish) 65 211.8811 388 987.950 1.28019 0.08307834 



WO 02/064820 



PCT/AU02/00122 



- 41 - 

Table 11 

Characteristics of the Cattle Sample 

Total: 166 87 

5 high instron compression 87 39 

low instron compression 89 38 



Breeds : 

Angus 

10 Belmont Red 

Brahman 
Hereford 
Santa Gertrudis 
Shorthorn 

15 

Regions : 

Pasture South 
Pasture North 
Grain South 
20 Grain North 



25 
25 
33 
32 
26 
25 



47 
39 
43 
37 



12 
12 
18 
10 
15 
17 



20 
21 
20 
26 



25 



Markets : 

Korean 

Domestic 

Japanese 



81 
47 
38 



22 
45 
22 



30 



Cohorts : 25 14 

Median: steers per cohort 6 3 

bottom quartile: 3 1 

top quartile: 10 11 



Sires : 



35 



Median: steers per sire 
bottom quartile: 
top quartile: 



113 
1 
1 
2 



62 
1 
1 
2 
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Table 12 

Distribution of LOX genotypes in the breeds in the combined 
sample . 

5 Breed Genotype 

11 12 22 

Angus 12 16 5 

Belmont Red 19 14 2 

Brahman 20 21 4 

10 Hereford 1 7 27 

Santa Gertrudis 18 5 1 

Shorthorn 23 11 3 
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Table 13 

Associations between LOX genotypes (loxg) and STIC. 
A. First Sample 

5 

Analysis of Deviance Table 
Gaussian model 
10 Response: stic 

Terms added sequentially (first to last) 



Df Deviance Resid.Df Resid. Dev F Value Pr(F) 



15 



NULL 










144 


66.41782 










market 




2 


4. 


68851 


142 


61.72932 


10 


.11346 


0. 


.0001124 


age 




1 


0. 


33061 


141 


61.39871 


1 


.42631 


0, 


.2356143 


cohort 




23 


25 


.50396 


118 


35.89474 


4 


.78382 


0. 


.0000000 


region 




3 


2. 


05111 


115 


33.84363 


2 


.94960 


0, 


.0371506 


breed %in% 


region 


8 


3. 


73686 


107 


30.10677 


2 


.01517 


0, 


.0537829 


loxg %in% 


{ region/breed ) 2 0 


9. 


94057 


87 


20.16620 


2 


.14425 


0, 


.0081786 



B. Second Sample 



25 Analysis of Deviance Table 



Gaussian model 
Response: stic 

Terms added sequentially < first to last) 
3 0 Df Deviance Resid.Df Resid. Dev F Value Pr(F) 



NULL 








77 


36.73118 










market 


2 


18 


.56015 


75 


18.17103 


50. 


57448 


0 


.0000000 


age 


1 


1 


.22303 


74 


16.94800 


6. 


66525 


0 


.0131586 


region 


3 


0 


.49228 


71 


16.45571 


0. 


89428 


0 


.4514796 


breed %in% region 


9 


2 


.26136 


62 


14.19435 


1. 


36933 


0 


.2303598 


loxg %in% (region/breed) 17 


5 


.93715 


45 


8.25720 


1. 


90331 


0 


.0433708 
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NOTE: cohort could not bo fitted as it required a model with more 
terms than degrees of freedom. cohort was dropped since that allowed 
a maximum of other terms to be fitted. 

C. Combined Sample 

Analysis of Deviance Table 



10 



15 



Gaussian model 
Response: stic 

Terms added sequentially (first to last) 

Df Deviance Resid.Df Resid. Dev F Value Pr(F) 



NULL 

market 

age 

cohort 
region 

breed %in% region 



2 9.28932 
1 0.66404. 

23 28.46167 

3 0.97412 
9 2.08684 



20 



loxg %±n% (region/breed) 24 12.75507 



222 105.0131 

220 95.7238 14.63402 0.0000015 

219 95.0598 2.09221 0.1500079 

196 66.5981 3.89890 0.0000002 

193 65.6240 1.02306 0.3840448 

184 63.5372 0.73056 0.6804229 

160 50.7821 1.67448 0.0327546 
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Table 14 



Associations between LOX genotypes (loxg) and STADH . 



A. First Sample 



10 



15 



Analysis of Deviance Table 
Gaussian model 
Response: stadh 

Terms added sequentially (first to last) 

Df Deviance Resid.Df Resid. Dev F Value Pr(F) 



NULL 










139 


4 


.118160 










market 




2 


0. 


045611 


137 


4 


.072549 


1 


.63665 


0 


.2008490 


age 




1 


1. 


030876 


136 


3 


.041674 


73 


.98193 


0 


.0000000 


cohort 




23 


1. 


131428 


113 


1 


.910246 


3 


.53035 


0 


.0000129 


region 




3 


0. 


053546 


110 


1 


.856699 


1 


.28094 


0 


.2863452 


breed %±n% 


region 


8 


0. 


192061 


102 


1 


.664639 


1 


.72293 


0 


.1050769 


loxg %in% 


(region/breed) 19 


0, 


508104 


83 


1 


.156535 


1 


.91919 


0 


.0229812 



20 



B. Second Sample 

Analysis of Deviance Table 



25 



30 



35 



Gaussian model 
Response: stadh 

Terms added sequentially (first to last) 

Df Deviance Resid.Df Resid. Dev F Value Pr(F) 

NULL 
market 
age 

region 

breed %in% region 







76 


3 


.362345 










2 


0.560443 


74 


2 


.801903 


25 


.18352 


0. 


0000000513 


1 


0.486387 


73 


2 


.315516 


43 


.71164 


0. 


0000000425 


3 


0.390962 


70 


1 


.924554 


11 


.71193 


0. 


0000090738 


9 


0.413632 


61 




1.510922 


4 


.13035 


0. 


0006637134 



loxg %in% (region/breed) 

17 1.021326 



44 



0.489595 5.39922 0.0000032684 
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C. Combined Sample 



Analysis of Deviance Table 



5 Gaussian model 
Response: stadh 



Terms added sequentially (first to last) 

Df Deviance Resid.Df Resid. Dev F Value Pr(F) 



10 NULL 

market 
age 

cohort 
region 
15 breed %in% region 



2 0.268508 
1 1.252352 

23 2.111912 

3 0.065737 
9 0.325818 



loxg %in% (region/breed) 

23 0.784128 



216 7.480742 

214 7.212234 7.78709 0.0005990 

213 5.959882 72.63985 0.0000000 

190 3.847971 5.32594 0.0000000 

187 3.782233 1.27098 0.2863505 

178 3.456415 2.09982 0.0326124 

155 2.672287 1.97746 0.0079540 
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Table 15 

Estimated sizes of effects of genotype substitutions at LOX on 
instron compression and adhesion of the semitendinosus muscle. 



Response : s tic . resid 



10 



Grand-mean 
Loxg 



15 



0.027133 se 

111 112 

0.03958 -0.02321 

se 0.07654 0.07897 



Response: sticbox3 . resid 

Grand-mean -0.012888 se 

Loxg Hi H2 

0.033317 -0.058769 
se 0.055916 0.058431 



0.053648 
122 
0.06503 
0.11751 



0.039137 
122 

-0.013213 
0.085115 



20 



Response : sticf ull . resid 

Grand-mean -0.012888 se 0.039137 

loxg in 112 122 

0.033317 -0.058769 -0.013213 

se 0.055916 0.058431 0.085115 



Response : sticadh . resid 
Grand-mean 0.002521 
25 Loxg in 

-0.001426 -0. 002304 
se 0.014860 0.015390 



se 
112 



0 .010321 
122 
0.011292 
0.022384 



Response : sticbox3adh . resid 
3 0 Grand-mean 0.004308 
loxg in 



se 
112 



-0.003154 -0.007446 
se 0.018594 0.018891 



0.012879 
122 
0.023524 
0 .028111 



35 
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Response: sticf \illadh. resid 

Grand-mean 0.0037399 se 0.0092835 

loxg Hi H2 122 

-0.009303 -0.005244 0.025767 

5 se 0.013379 0.013763 0.020180 
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Table 16 



Information on the data sets 



Effect 


Class 


LDIC 


LDPF 


Combined 
Extreme 


Non- 
extreme 


Combined 


Total 




136 


131 


398 


543 


916 


Sires 




96 


96 


171 


61 


227 


Cohorts 




24 


24 


26 


10 


30 


Breeds 


Angus 


19 


22 


62 


134 


196 




Belmont Red 


23 


25 


66 


140 


200 




Brahman 


24 


27 


76 


73 


142 




Hereford 


27 


14 


67 


1 


68 




Santa Gertrudis 


25 


24 


73 


195 


257 




Shorthorn 


18 


19 


53 


0 


53 


Regions 


Pasture South 


27 


27 


92 


93 


185 




Pasture North 


24 


25 


78 


144 


212 




Grain South 


54 


41 


125 


168 


291 




Grain North 


31 


38 


102 


13 8 i 


228 


Markets 


Domestic 


47 


51 


130 


240 


361 




Korean 


61 


58 


189 


201 


376 




Japaness 


28 


22 


78 


102 


179 
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Table 17 



Meat quality traits tested for LOX gene marker 



Code 


Trait 


LD Fat% 


Tn^ramncnil ar Pa f- nprronhAna f CnvfiTrl of* Mof hn/^\ 


LD_IC 


T.oTi rr i q q i mil q f^m^qi Tn c ^ Trtm *~\t\\*"\ >* q g a ^ r*r* 
uviiy x o s xmus uuxax xuo umiu uUiitproSaioiJ 


LD IY 


iiOflyissiiuus uorsi lnitiai yieia Kins only/ 


LD LOSS 


uunyisbimus uoxm luujvi. ny losses 


LD PF 


ijouyissiinu.3 aoisx roajt r ore © -~m\is c use ixn aiso 


LD PF-IY 


iionyis s iiuus aorsi reaA r orce — initial yield (Ntn.) 




LonQTissimus dorsi a* colour 


Li/ JJ 


Longissinvus dorsi b* colour 


t n 1 


Longiss imus dorsi L* colour 


T n i-»w 
■Li -U pit 


Longissimus dorsi. ultimate pH 


w x xv i? a u 'o 


Intramuscular Fat percentage (NIR method) 


ST AdhRS 




ST_IC 


Semitendinosus Instrom compression 


ST_IY 


Semitendinosus initial yield (Nth kills only) 


ST_LOSS 


Semitendinosus cooking loss% 


ST_PF 


Semitendinosus Peak Force 


ST_PF-IY 


Semitendinosus Peak Force - initial yield (Nth) 


ST_a 


Semitendinosus a* colour 


ST_b 


Semitendinosus b* colour 


ST_1 


Semitendinosus L* colour 


ST_pH 


Semitendinosus ultimate pH 


TenderQ 


Tenderness Quality as measured by PF (x 100) 



5 
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Table 18 

Distribution of LOX genotypes in the breeds in the three datasets 



Breed 


Extreme 


Non- extreme C 


Combined 




11 


12 


22 


11 


12 


22 


11 


12 


22 


Angus 


21 


33 


7 


38 


75 


21 


59 


109 


28 


Brahman 


35 


35 


6 


19 


40 


14 


52 


71 


19 


Belmont Red 


31 


23 


12 


59 


60 


21 


87 


83 


30 


Hereford 


4 


17 


46 


0 


1 


0 


4 


18 


46 


Santa Gertrudis 


42 


29 


2 


120 


62 


13 


156 


86 


15 


Shorthorn 


34 


18 


1 


0 


0 


0 


34 


18 


1 


Total 


167 


156 


74 


236 


238 


69 


392 


385 


139 ' 



5 
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Table 19 



10 



15 



Dependent Variable: STIC 



Source 



<.0001 



Model 

1 

Error 

Corrected Total 



297 



The GLM Procedure 



Sum of 
Squares 



90.3108958 



97 11.5922526 
394 101.9031484 



Mean Square F Value 



0.3040771 



0.1195078 



2.54 



20 



25 



R-Square 
0.886242 



Source 



Coeff Var 
16.41081 



Root MSE 
0.345699 



STIC Mean 
2.106532 



DF 



Type III SS 



Mean Square 



F Value 



Pr 



contemp ( Bcode ) 



30 



35 



<.0001 



Fingp (Bcode) 



0.2004 



Stim 



0.1199 



lox( Bcode) 



0.0145 



0.0243 



269 
3 
1 
12 
1 



72.05453149 
0.56463473 
0.29426415 
3.23498136 
0.62547588 



0 .26786071 
0.18821158 
0.29426415 
0.26958178 
0.62547588 



2.24 
1.57 
2.46 
2.26 
5.23 



40 



Contrast 



Contrast SS 



Mean Square 



F Value 



Additive Test 11-22 



45 



50- 



55 



60 



65 



70 



0.0636 



11-12 



<.0001 



12-22 



0.2283 



Dominance Test 



0.0052 



Recessive Test 



0.7770. 



Over Dominance Test 



0.0004 



Parameter 

Additive Test 11-22 

11- 12 

12- 22 

Dominance Test 
Recessive Test 
Over Dominance Test 



LDL 



Dependent Variable: LD1 



0.42074192 
2.17031475 
0.17567132 
0.97536352 
0.00964320 
1.5877 6355 



Estimate 

-0.54698783 
-0.93226650 
0.38527867 
1.31754517 
0.16170916 
-1.47925433 



0.42074192 
2.17031475 
0. 17567132 
0.97536352 
0.00964320 
1.S8776355 



Standard 
Error 

0.29151964 
0.21876446 
0.31777709 
0.46119045 
0.56927504 
0.40583359 



3.52 
18. 16 
1.47 
8. 16 
0.08 
13.29 



t Value 

-1.88 
-4.26 
1.21 
2.86 
0.28 
-3.64 



Pr > |t| 

0.0636 
<.0001 
0.2283 
0.0052 
0.7770 
0.0004 



The GLM Procedure 
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Table 19 (Continued) 



10 



15 



20 



25 



30 



35 



40 



45 



50 



55 



Source 



Model 



< . 0001 



Error 

Corrected Total 



R-Square 
0.857745 



Source 



contemp (Bcoda) 



0.0278 



Fingp(Bcode) 



0.3161 



Stira 



0.8678 



lox(Bcode) 



0.0427 



wt 



0.3590 



Contrast 



Additive Test 



0.3625 



11 vs 12 



0.5091 



12 vs 22 



0.7024 



Dominance Test 



0.9600 



Recessive Test 



0.4966 



Over Dominance Test 



0.3129 



295 

96 
391 



Sum of 
Squares 



4070.356476 

675.059834 
4745.416310 



Mean Square F Value Pr > 
13.797819 1.96 
7.031873 



Coeff Var 
6.920726 



Root MSE 
2.651768 



LD1 Mean 
38.31633 



DF 

267 
3 
1 
12 
1 

DF 

1 
1 
1 
1 
1 
1 



Type III SS 

2628.383772 
25.195055 
0.195851 
160.954945 
5.974613 

Contrast SS 

5 . 88718974 
3.08780445 
1.03256209 
0.01778689 
3.27515251 
7.23682011 



Mean Square 

9.844134 
8.398352 
0. 195851 
13.412912 
5.974613 

Mean Square 

5.88718974 
3.08780445 
1.03256209 
0.01778689 
3.27515251 
7.23682011 

Standard 



F Value 

1.40 
1.19 
0. 03 
1.91 
0.85 

F Value 

0.84 
0.44 
0.15 
0.00 
0.47 
1.03 



Pr > 



Pr > 



Parameter 




Estimate 




Error 


t Value 


Pr > |t| 


Additive Test 11-22 


-2 


. 04608559 


2 . 


.23617248 


-0.91 


0.3625 


11-12 


-1, 


. 11200497 


1. 


.67809812 


-0.66 


0 .5091 


12-22 


-0, 


. 93408062 


2 , 


.43759636 


-0.38 


0 . 7024 


Dom inane e Te s t 


0, 


.17792435 


3 


.53769859 


0.05 


0.9600 


Recessive Test 


2, 


. 98016622 


4 . 


.36676923 


0.68 


0. 4966 


Over Dominance Test 


-3, 


.15809056 


3 . 


.11305081 


-1.01 


0.3129 
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Table 20 (continued) 

Single Trait Logistic Regression 



10 



15 



20 



25. 



30 



35 



The LOGISTIC Procedure 
Model Information 



Data Set 

Response Variable 
Number of Response Levels 
Number of Observations 
Link Function 
Optimization Technique 



WORK . EXTRERES 

lox 

3 

389 
Logit 

Fisher's scoring 



Response Profile 



Ordered 



Total 



Value 


lox 


Frequency 


1 


11 


165 


2 


12 


153 


3 


22 


71 



NOTE: 8 observations were deleted due to missing values for the response or explanatory- 
variables. 

The LOGISTIC Procedure 
Testing Global Null Hypothesis: BETA=0 
Test Chi-Square DF Pr > ChiSq 



Likelihood Ratio 

Score 

Wald 



5.0720 
5.0502 
5.0083 



0.0243 
0.0246 
0.0252 



40 



45 



Parameter 

Intercept 

Intercept2 

sticpred 



Analysis of Maximum Likelihood Estimates 

Chi-Square Pr > ChiSq 



DF 



Estimate 



Standard 
Error 



0.9245 
2.7482 
-0.5827 



0.5608 
0.5777 
0.2604 



2.7174 
22.6278 
5.0083 



0.0993 
<.0001 
0.0252 



50 



Odds Ratio Estimates 



Effect 
sticpred 



, Point 
Estimate 



0.558 



95% Wald 
Confidence Limits 



0.335 



0.930 
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Table 20 (continued) 



Multi-trait Logistic Regression 



10 



15 



20 



25 



Test 



The LOGISTIC Procedure 
Testing Global Null Hypothesis: BETA=0 

Chi-Square DF Pr > ChiSq 



Likelihood Ratio 

Score 

Wald 



14.7234 
14.5859 
13.9424 



0.0021 
0.0022 
0.0030 



Analysis of Maximum Likelihood Estimates 











Standard 


Wald 




Parameter 




DF 


Estimate 


Error 


Chi-Square 


Pr > ChiSq 


Intercept 


11 


1 


-0.9210 


0.8702 


1.1202 


0.2899 


Intercept 


12 


1 


0.9467 


0.8711 


1.1813 


0.2771 


sticpred 




1 


-0.8G81 


0.3191 


7.4014 


0.0065 


stadhpred 




1 


-0.4025 


0.7724 


0.2715 


0.6023 


stpf pred 




1 


0.5689 


0.1936 


8.6313 


0.0033 



Odds Ratio Estimates 

Point 95% Wald 

Effect Estimate Confidence Limits 

sticpred 0.420 0.225 0.785 

stadhpred 0.669 0.147 3.039 

2.581 



30 



35 stpfpred 1.766 1.208 
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Table 21 



10 



15 



The Mixed Procedure 
Model Information 



Data Set 

Dependent Variable 
Covariance Structure 
Estimation Method 
Residual Variance Method 
Fixed Effects SE Method 
Degrees of Freedom Method 



WORK.CRC1 
STIC 

Variance Components 

REML 

Profile 

Model-Based 

Containment 



Covariance Parameter Estimates 



20 



25 



Cov Parm 

SirelD(Bcode) 
contemp {Bcode) 
Residual 



Estimate 

0.002544 
0.008813 
0.07009 



Standard 
Error 

0.002147 
0.003896 
0.005174 



Z 

Value 

1.18 
2.26 
13 .55 



Pr 2 

0.1181 
0.0118 
<.0001 



Type 


3 Tests 


of Fixed 


Effects 






Num 


Den 






Effect 


DF 


DF 


F Value 


Pr > F 


Fingp(Bcode) 


9 


317 


17.50 


<.0001 


Stim 


1 


317 


1 .75 


0. 1870 


lox (Bcode) 


8 


317 


2.3 6 


0.0176 


wt 


1 


317 


6.45 


0.0115 



Estimates 







Standard 










Label 


Estimate 


Error 


DF 


t Value 


Pr 


> Itl 


Additive Test 11-22 


0.05547 


0.07813 


317 


0.71 


0 


.4782 


11-12 


0.1044 


0.05660 


317 


1.84 


0 


.0661 


12-22 


-0.04890 


0.07013 


317 


-0.70 


0 


.4861 


Dominance Test 


-0.1533 


0.1007 


317 


-1.52 


0 


. 1289 


Recessive Test 


-0.00657 


0.1373 


317 


-0.05 


0 


.9619 


OverDominance Test 


0.1598 


0.1170 


317 


1.37 


0 


.1730 



50 

Contrasts 





Num 


Den 






Label 


DF 


DF 


F Value 


Pr > F 


Additive Test 11-22 




317 


0.50 


0^4782 


11-12 




317 


3.40 


0.0661 


12-22 




317 


0.49 


0.4861 


Dominance Test 




317 


2.32 


0.1289 


Recessive Test 




317 


0.00 


0.9619 


OverDominance Test 




317 


1.87 


0.1730 



lnstramuscuar Fat 

65 

The Mixed Procedure 

Model Information 

Data Set W0RK.CRC1 
Dependent Variable Fat 

Covariance Structure Variance Components 
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Table 21 (continued) 



Estimation Method REML 

Residual Variance Method Profile 

Fixed Effects SE Method Model-Based 

Degrees of Freedom Method Containment 



10 



15 



Covariance Parameter Estimates 



Cov Parm 

SirelD(Bcode) 
contemp (Bcode) 
Residual 



Estimate 

0 . 04848 
0.3111 
1 .2479 



Standard 
Error 

0.04510 
0.1060 
0.09791 



2 

Value 

1 .07 
2.93 
12.75 



Pr Z 

0.1412 
0.0017 
<.0001 



20 



Effect 



Type 3 Tests of Fixed Effects 



Num 
DF 



Den 
DF 



F Value 



Pr > F 



25 



Fingp (Bcode) 
Stim 

lox (Bcode) 
wt 



306 
306 
306 
306 



7.08 
0.31 
2.00 
74.54 



<.00C1 
0.5790 
0.0461 
<.0001 



30 



Estimates 



35 



40 



Standard 

Label Estimate Error DF t Value 

Additive Test 11-22 0.3765 0.3427 306 1.10 

11- 12 -0.1207 0.2505 306 -0.48 

12- 22 0.4972 0.3029 306 1.64 
Dominance Test 0.6179 0.4377 306 1.41 
Recessive Test -0.8737 0.5964 306 -1.47 
OverDominance Test 0.2558 0.5184 306 0.49 



Pr 



0.2728 
0.6304 
0.1018 
0.1591 
0.1439 
0.6220 



45 



Label 



Contrasts 



Num 
DF 



Den 
DF 



F Value 



Pr 



Additive Test 11-22 


1 306 


1 


21 


0 


2728 


11-12 


1 306 


0 


23 


0 


6304 


12-22 


1 306 


2 


69 


0 


1018 


Dominance Test 


1 306 


1 


99 


0 


1591 


Recessive Test 


1 306 


2 


15 


0 


1439 


OverDominance Test 


1 306 


0 


24 


0 


6220 



55 

LDPH 

The Mixed Procedure 



60 



65 



Model Information 



Data Set 

Dependent Variable 
Covariance Structure 
Estimation Method 
Residual Variance Method 
Fixed Effects SE Method 



WORK.CRC1 
LDpH 

Variance Components 

REML 

Profile 

Model-Based 



Degrees of Freedom Method Containment 



7 0 Covariance Parameter Estimates 

Standard 



WO 02/064820 PCT/AU02/00122 

- 58 - 



Table 21 (continued) 



Cov Parm Estimate Error Value Pr Z 

5 contempt Bcode) 0.002771 0.000702 3.95 <.0O01 

Residual 0.007878 0.000572 13.77 <.0001 



Type 


3 Tests 


of Fixed 


Effects 






Num 


Den 






Effect 


DF 


DF 


F Value 


Pr > F 


Fingp(Bcode) 


9 


363 


3 .44 


0.0004 


Stim 


1 


363 


1.14 


0.2854 


lox (Bcode) 


8 


363 


3.01 


0.0027 


wt 


1 


363 


21.61 


<.O001 



2 0 Estimates 

Standard 



Label 


Estimate 


Error 


DF 


t 


Value 


Pr > |c| 


Additive Test 11-22 


-0. 1089 


0.02648 


363 




-4.11 


< 


.0001 


11-12 


-0.01502 


0.01922 


363 




-0.78 


0 


.4351 


12-22 


-0.09393 


0.02376 


• 363 




-3.95 


< 


.0001 


Dominance Test 


-0.07891 


0.03416 


363 




-2.31 


0 


.0215 


Recessive Test 


0.2029 


0.04650 


363 




4.36 


< 


.0001 


OverDominance Test 


. -0.1240 


0.03971 


363 




-3.12 


0 


0019 




The Mixed Procedure 
















Contrasts 
















Num Den 












Label 




DF DF 


F Value 


Pr > F 






Additive Test 


11-22 


1 363 


16 


93 


<.0001 






11-12 




1 363 


0 


61 


0.4351 






12-22 




1 363 


15 


63 


<.0001 






Dominance Test 




1 363 


5 


34 


0.0215 






Recessive Test 




1 363 


19 


04 


<.0001 






OverDominance 


Test 


1 363 


9 


75 


0.0019 
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Table 22 
STIC 

The Mixed Procedure 
Covariance Parameter Estimates 



10 







Standard 


Z 






Cov Parm 


Estimate 


Error 


Value 




Pr Z 


SirelD(Bcode) 


0.006141 


0.004815 


1.28 


0 


1011 


contemp (Bcode) 


0.05433 


0.01003 


5.42 


< 


0001 


Residual 


0.09190 


0.006489 


14.16 


< 


0001 



20 Effect 



Type 3 Tests of Fixed Effects 

Num Den 
DF DF F Value Pr > F 



Fingp (Bcode) 


11 


369 


9 


24 


< 


0001 


Stim 


2 


369 


4 


06 


0 


0181 


lox (Bcode) 


12 


369 


1 


90 


0 


0336 


wt 


1 


369 


9 


93 


0 


0018 



Estimates 



30 



35 



40 







Standard 




Label 


Estimate 


Error 


DF 


Additive Test 11-22 


0.06327 


0.08142 


369 


11-12 


0.07985 


0.05752 


369 


12-22 


-0.01658 


0.07440 


369 


Dominance Test 


-0.09643 


0.1052 


369 


Recessive Test " 


-0.04669 


0.1450 


369 


OverDominance Test 


0.1431 


0.1198 


369 



t Value Pr > |t[ 



0.78 



39 
22 
92 



-0.32 
1.20 



4376 
1659 
8238 
3598 
7476 
2328 



Contrasts 



45 



50 



Label 



Num 
DF 



Den 
DF 



F Value 



Pr > F 



Additive Test 11-22 


1 369 


0 


60 


0 


4376 


11-12 


1 369 


1 


93 


0 


1659 


12-22 


1 369 


0 


05 


0 


8238 


Dominance Test 


1 369 


0 


84 


0 


3598 


Recessive Test 


1 369 


0 


10 


0 


7476 


OverDominance Test 


1 369 


1 


43 


0 


2328 



STL 

55 Covariance Parameter Estimates 

Standard 2 





Cov Parm 


Estimate 


Error 


Value 




Pr Z 


60 


SirelD(Bcode) 


0.4389 


0.3429 


1.28 


0 


1003 




contemp (Bcode) 


3.3608 


0.7375 


4.56 


< 


0001 




Residual 


10.6090 


0.6789 


15.63 


< 


0001 



Type 


3 Tests 


of Fixed 


Effects 






Num 


Den 






Effect 


DF 


DF 


F Value 


Pr > F 


Fingp (Bcode) 


.11 


368 


16.40 


<.0001 


Stim 


2 


368 


4.68 


0.0098 


lox (Bcode) 


12 


368 


2.18 


0.0124 
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Table 22 (Continued) 



- 60 - 



wt 



368 



8.64 0.0035 



Estimates 



Standard 



10 



15 



20 



25 



30 



Label 


Estimate 




Error 


DF 


t 


Value 


Pr > |c| 


Additive Test 11-22 


-0.7500 




0.8582 


368 




-0 . 87 


0. 


.3827 


11-12 


-0.6732 




0.5998 


368 




-1 . 12 


0 


.2625 


12-22 


-0.07683 




0.7815 


368 




-0.10 


0. 


.9217 


Dominance Test 


0.5963 




1.0975 


368 




0.54 


. 0. 


.5872 


Recessive Test 


0.8268 




1.5279 


368 




0.54 


0 


5887 


OverDomi nance Test 


-1.4232 




1.2576 


368 




-1.13 


0. 


.2585 






Contrasts 
















Num 


Den 












Label 




DF 


DF 


F Value 


Pr > F 






Additive Test 


11-22 


1 


368 


0 


.76 


0.3827 






11-12 




1 


368 


1 . 


.26 


0.2625 






12-22 




1 


368 


0 


.01 


0.9217 






Dominance Test 


1 


368 


0. 


.30 


0.5872 






Recessive Test 




1 


368 


0. 


.29 


0.5887 






OverDomi nance 


Test 


] 


368 


1 . 


.28 


0.2585 
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Table 23 



The Mixed Procedure 



Covariance Parameter Estimates 



10 



15 



20 



25 



30 



35 







Standard 




Z 










Cov Parm 


Estimate 


Error 


Value 




Pr Z 






SirelD(Bcode) 


0.005452 


0.004008 


1 


.36 


0. 


.0869 






contemp (Bcode) 


0.03915 


0.008774 


4 


.46 


< . 


.0001 






Residual 


0.09105 


0.006426 


14 


. 17 


< . 


.0001 








Type 3 Tests 


of Fixed Effects 














Num 


Den 














Effect 


DF 


DF F 


Value 


Pi 


> F 








Fingp(Bcode) 10 


368 


10.46 


< . 


0001 








Stim 


2 


368 


5.89 


0. 


0030 








lox (Bcode) 


10 


368 


2 .28 


0. 


0132 








wt 


1 


368 


7.14 


0. 


0079 










Estimates 


















Standard 














Label 


Estimate 


Error 


DF 


t 


Value 


Pr > |t[ 


Additive Test 11-22 


0.06263 


0.07965 


368 




0.79 




0. 


.4322 


11-12 


0.09181 


0.05618 


368 




1.63 




0. 


1031 


12-22 


-0.02918 


0.07277 


368 




-0.40 




0. 


.6887 


Dominance Test 


-0.1210 


0.1028 


366 




-1.18 




0. 


2398 


Recessive Test 


-0.03345 


0.1419 


368 




-0.24 




0. 


8137 


OverDominance Test 


0.1544 


0.1171 


368 




1.32 




0. 


1879 



40 



45 



Contrasts 





Num 


Den 






Label 


DF 


DF 


F Value 


Pr > F 


Additive Test 11-22 


1 


368 


0.62 


0.4322 


11-12 


1 


368 


2.67 


0. 1031 


12-22 


1 


368 


0.16 


0. 6887 


Dominance Test 


1 


368 


1.39 


0.2398 


Recessive Test 


1 


368 


0.06 


0.8137 


OverDominance Test 


1 


368 


1.74 


0.1879 
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Claims : 

1. A method for assessing the tenderness of meat from an 
animal, comprising the step of testing the animal for the 

5 presence or absence of a genetic marker selected from the 
group consisting of: 

(1) an allele of the gene encoding calpastatin 
(CAST) associated with peak-force variation or 
genetic variation located other than in the CAST 

10 gene which shows allelic association with the 

CAST allele; and 

(2) an allele of the gene encoding lysyl oxidase 
(LOX) associated with variation in instron 
compression of the semitendinosis muscle or 

15 genetic variation located other than in the LOX 

gene which shows allelic association with the 
LOX allele. 

2. A method as claimed in claim 1 wherein the genetic 
2 0 marker is located in the 3 ' UTR of the CAST gene . 

3. A method as claimed in claim 2 wherein the allele is 
CAST3 D/E allele 1. 

25 4. A method as claimed in claim 2 wherein the allele is 
CAST 3 D/E allele 2 or allele 3. 



5. A method as claimed in 
marker is located in the 5' 

30 

6* A method as claimed in 
CAST 5 allele 3. 

7. A method as claimed in 
3 5 CAST 5 allele 7 or allele 9. 

8. A method as claimed in 



claim 1 wherein the genetic 
UTR of the CAST gene. 

claim 5 wherein the allele is 
claim 1 wherein the allele is 
claim 1 wherein the genetic 
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marker is located in the genomic region embracing the 
coding sequence of the CAST gene. 

9. A method as claimed in claim 1 wherein the genetic 
5 marker is a genetic variation located other than in the 
CAST gene which shows allelic association with either or 
both of the CAST 3 D/E polymorphism and the CAST 5 
polymorphism . 

10 10. A method as claimed in any one of claims 1 to 9, 

further comprising the step of testing for the presence or 
absence of one or more additional genetic markers 
associated with peak- force variation. 

15 11. A method as claimed in claim 10 further comprising 

the step of testing for the presence of CAST5 allele 7 or 
allele 9 and/or the absence of CASTS allele 3 once the 
presence of CAST D/E allele 2 has been established. 

20 12. A method as claimed in claim 1 wherein the genetic 
marker is allele 1 of the LOX polymorphism. 

13. A method as claimed in claim 1 wherein the genetic 
marker is allele 2 of the LOX polymorphism. 

25 

14. A genetic marker for meat tenderness in an animal 
which is a polymorphic form of the CAST gene, being the 
CAST 3 D/E polymorphism or the LOX polymorphism. 

3 0 15. A genetic marker as claimed in claim 14 wherein the 

variable portion of the CAST gene comprises the nucleotide 
sequence set forth in SEQ ID N0:1. 

16. A genetic marker as claimed in claim 14 wherein the 
3 5 variable portion of the CAST gene comprises the nucleotide 
sequence set forth in SEQ ID NO: 2. 
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17. A genetic marker as claimed in claim 14 wherein the 
variable portion of the CAST gene comprises the nucleotide 
sequence set forth in SEQ ID NO: 3. 

5 

18. An isolated DNA molecule comprising the nucleotide 
sequence set forth in SEQ ID NO:l. 

19. An isolated DNA molecule comprising the nucleotide 
10 sequence set forth in SEQ ID NO: 2. 

20. An isolated DNA molecule comprising the nucleotide 
sequence set forth in SEQ ID NO: 3. 

15 21. An isolated DNA molecule consisting of the nucleotide 
sequence set forth in SEQ ID NO:l. 

22. An isolated DNA molecule consisting of the nucleotide 
sequence set forth in SEQ ID NO: 2. 

20 

23. An isolated DNA molecule consisting of the nucleotide 
sequence set forth in SEQ ID NO: 3. 

24. A method for selecting an animal likely to yield meat 
25 of improved tenderness, comprising the steps of: 

(1) testing the animal for the presence of an allele 
of the gene encoding calpastatin (CAST) 
associated with low peak-force or genetic 
variation located other than in the CAST gene 

3 0 which shows allelic association with the CAST 

allele, and/or for the presence of an allele of 
the LOX gene associated with low instron 
compression of the semitendinosis muscle or 
genetic variation located other than in the LOX 

3 5 gene which shows allelic association with the 

LOX allele; and 

(2) selecting animals which have the CAST and/or LOX 
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allele and/or genetic variation in allelic 
association therewith. 

25. A method as claimed in claim 24 wherein the allele 
5 tested for is CAST3 D/E allele 2 and animals which are 

homozygous for this allele are selected, 

26. A method as claimed in claim 25 further comprising 
the step of testing for the presence or absence of one or 

10 more additional alleles of the gene encoding calpastatin 
associated with low peak-force. 

27. A method as claimed in claim 26 wherein the presence 
of CASTS allele 7 or allele 9 is tested for, and animals 

15 having either of these alleles in addition to CAST 3 D/E 
allele 2 are selected. 



28. A method as claimed in claim 26 wherein the presence 
of CASTS allele 3 is tested for, and animals having this 

20 allele are rejected despite the presence of CAST3 D/E 
allele 2. 

29. A method as claimed in any one of claims 24 to 27, 
further comprising using the selected animal for breeding. 

25 

30. An oligonucleotide probe for amplification of a 
genetic marker associated with peak-force variation, said 
genetic marker being either an allele of the gene encoding 
calpastatin (CAST) or genetic variation located other than 

3 0 in the CAST gene which shows allelic association with said 
allele. 

31. An oligonucleotide probe as claimed in claim 30 
selected from the group consisting of: 

3 5 castd 5' cat ttg gaa aac gat gcc tea c 3' 

caste 5' tct acg att age age tea aga gga g 3' 
CAST5U1 5 ' - GTAAAGCCGC AC AAAACACACCC AGG- 3 ' 
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CAST5D1 5 ' -GTTTCTGGACCCTCTGGATGAGGAAGCGG- 3 ' 

32. An oligonucleotide probe for amplification of a 
genetic marker associated with variation in instron 

5 compression of the semitendinosis muscle, the genetic 

marker being either an allele of the gene encoding lysyl 
oxidase (LOX) or genetic variation located other than in 
the LOX gene which shows allelic association with said 
allele. 

10 

33. An oligonucleotide probe as claimed in claim 31 
selected from the group consisting of: 

LOX K5 : 5 ' tat cac tga tgt caa acc tg 3 ' 
LOX K6: 5' act cag gca cca aat age tg 3' 

15 

34 . A kit for use in assessing the tenderness of meat 
from an animal and/or selecting an animal likely to yield 
meat of improved tenderness, comprising oligonucleotide 
probes for amplification of at least one genetic marker 

20 associated with meat tenderness, said genetic marker being 
either an allele of the gene encoding calpastatin (CAST) 
or genetic variation located other than in the CAST gene 
which shows allelic association with said allele, or an 
allele of the LOX gene associated with low instron 

2 5 compression of the semitendinosis muscle or genetic 

variation located other than in the LOX gene which shows 
allelic association with the LOX allele, and means for 
amplifying DNA. 

3 0 35. A kit as claimed in claim 34 wherein the 

oligonucletide probes are selected from the group 

consisting of: 

castd 5' cat ttg gaa aac gat gec tea c 3' 

caste 5' tct acg att age age tea aga gga g 3 

3 5 CAST5U1 5 ' - GTAAAGCCGCACAAAAC ACACCCAGG- 3 ' 

CAST5D1 5 ' -GTTTCTGGACCCTCTGGATGAGGAAGCGG- 3 ' 
LOX K5: 5' tat cac tga tgt caa acc tg 3' 
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LOX K6: 5' act cag gca cca aat age tg 3' 

36. An animal when selected by the method as defined in 
any one of claims 24 to 27. 

5 

37. The progeny of an animal as defined in claim 36. 

38. Meat from an animal as defined in claim 36. 

10 39. Meat from the progeny of an animal as defined in 
claim 36. 

40. The use of an animal as defined in claim 36 in 
breeding . 



15 



41. The use of the progeny of an animal as defined in 
claim 3 6 in breeding. 
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170 

<210> 2 
<211> 170 
<212> DNA 
<213> Bos sp . 

<400> 2 
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